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Abstract 

We consider boundary conditions compatible with discrete holomorphicity for the dilute 
O(n) and loop models. In each model, for a general set of boundary plaquettes, multiple 
types of loops can appear. A generalisation of Smirnov's parafermionic observable is there- 
fore required in order to maintain the discrete holomorphicity property in the bulk. We show 
that there exist natural boundary conditions for this observable which are consistent with in- 
tegrability, that is to say that, by imposing certain boundary conditions, we obtain a set of 
linear equations whose solutions also satisfy the corresponding reflection equation. In both loop 
models, several new sets of integrable weights are found using this approach. 

1 Introduction 

The mathematically rigorous understanding of two-dimensional lattice models at criticality and 
their scaling limits has seen some significant recent breakthroughs, see e. g. [UEII3]. An impor- 
tant tool in this undertaking is the discretely holomorphic parafermionic observable introduced 
by Smirnov PQ. The premise for the existence of such observables is as follows: at the continuum 
level, when the lattice spacing vanishes, many lattice models are believed to be conformally in- 
variant and admit holomorphic observables. One therefore expects that there exist observables 
in the finite lattice models that satisfy a discrete form of holomorphicity. 

Discretely holomorphic observables have the property that their discrete contour integral 
around any closed path vanishes. A similar property has been discussed by Bernard and Felder 
in the context of non-local conserved currents in lattice quantum field theories [3]. However, 
the true power of discretely holomorphic observables was shown by Smirnov in proving the 
conformal invariance of the Ising model [5j. They have since been used by Duminil-Copin and 
Smirnov to obtain a rigorous proof of the connective contstant of self-avoiding walks on the 
honeycomb lattice [6], a value predicted 30 years earlier by Nienhuis [7J. Other parafermionic 
observables and connections with SLE are discussed in [HI ED] • 
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There is a remarkable but still not fully understood connection between discrete holomor- 
phicity and Yang-Baxter integrability: by choosing a contour of integration surrounding a single 
plaquette of the lattice, one readily obtains the integrable Boltzmann weights of the model. This 
connection with discrete holomorphicity has been established by Rajabpour and Cardy [IT] in 
the Zn model, and by Ikhlef and Cardy [121 [13] m the Potts model, the dense and dilute 0{n) 
models and the loop model. For a discussion of discrete holomorphicity in the six- and 
eight- vertex models, see [13] . In the case of the Z^ model, Alam and Batchelor [15J have recently 
shown that considering a discrete contour integral around two or three adjacent plaquettes leads 
to the inversion and star-triangle equations, respectively. We remark here that the work [I] of 
Bernard and Felder mentioned above discusses integrability in the context of quantum groups 
and the related R-matrices, and that Nienhuis [T6l [TT] derived integrable weights for the 0{n) 
loop models from a criticality argument. 

An obvious question is whether there exist boundary conditions for discretely holomorphic 
observables that give rise to integrable boundary weights, i.e. solutions to the reflection equa- 
tion [18] , a boundary version of the Yang-Baxter equation. Systematic studies of such solutions 
for interaction-round-a-face models and loop models are found in [19J and [20| [21] . 

In [22], Beaton et al used a variant of discrete holomorphicity at the boundary to rigor- 
ously prove the critical temperature of the adsorption transition of self-avoiding walks on the 
honeycomb lattice, a value that had been obtained by Batchelor and Yung [23J by solving the 
corresponding reflection equation. Following the same reasoning as in [22], Beaton has also con- 
sidered [21] a rotated honeycomb lattice with a boundary. Ikhlef imposed [22] certain boundary 
conditions on a discretely holomorphic observable in the 0(n) and Zjy models. In the 0{n) 
case, he recovered the integrable diagonal weights, and in the Zjsi case he found new integrable 
weights. 

The aim of this paper is to further the understanding of the connection between boundary 
integrability and boundary conditions of discretely holomorphic observables in the case of more 
general non-diagonal boundary weights, where loops attached to the boundary may acquire 
different weights compared to loops in the bulk, such as in j2HH2Z]- In the process, we slightly 
generalise the observables defined in the 0{n) and loop models and show they do indeed 
satisfy natural boundary conditions compatible with integrability. These boundary conditions 
give rise to a set of linear equations whose solutions satisfy the corresponding reflection equa- 
tions. In this way, we not only recover the known solutions to the reflection equations, but also 
obtain new ones. 

The structure of this paper is as follows. In Section [2j we review the dilute 0(n) model and 
the definition of the associated discretely holomorphic parafermionic observable. In Section [3j 
we extend this definition to ensure that this observable remains discretely holomorphic in the 
presence of non-trivial boundaries. In Section [4j we introduce two boundary conditions that the 
observable can satisfy, and use these to obtain the integrable boundary weights of the model. In 
Section^, we repeat this procedure for the loop model. Section ^ contains some concluding 
remarks, while a discussion of the reflection equations and a further generalisation of the dilute 
O(n) model are deferred to Appendix |A| 



2 The dilute O(n) loop model 

2.1 Description of the model 

The dilute 0(n) model is a lattice model of closed non-intersecting, non-oriented loops. A given 
configuration of the model consists of a choice of tiling of the domain ft (see below) by the bulk 
plaquettes shown in Fig. [I] and the boundary plaquettes shown in Fig. [2] 
Four types of loops can occur in the model. 
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Figure 1: The bulk plaquettes of the dilute 0(n) model. The Boltzmann weights are 
t,u\,U2,w\,W2 and v, respectively. 




Figure 2: The boundary plaquettes of the dilute 0(n) model. The Boltzmann weights 
are /?i,/?2 and /J3, respectively. 



• Bulk loops that lie entirely in the bulk. These loops carry the weight n. 

• Bulk loops that pass through the boundary via 02 plaquettes. These loops also carry the 
weight n. 

• Boundary loops that are attached to the boundary at a single /3s boundary plaquette, or 
at two different (3% plaquettes in such a way that the top part of the loop lies in the top 
half of the topmost of the two boundary plaquettes. These loops receive the weight n\. 

• Boundary loops that are attached to the boundary at two different /3s boundary plaquettes 
such that the top half of the loop lies in the bottom half of the topmost of the two boundary 
plaquettes. These loops receive the weight 77,2. 

Let us also define the domain fi: starting from a regular square lattice, we skew adjacent 
rows in an opposite sense and scale the domain so all the faces become rhombi with acute angle 
a. There is also a right-hand boundary consisting of isosceles triangles. A small portion of 
is shown in Fig. [3j For a = |, the domain becomes the regular square lattice with a single 
straight-edged boundary. 

The total weight P(G) of a configuration G is the product of three contributions: 

P{G) = ^#[^-P lac l uettes ] ^#[i;-plaquettes]^ ^#[/3i-plaquettes]^#[/3 2 -plaquettes]^#[/3 3 -plaquettes] ^ 

x ( n # [ n - lo °P s ] n f [ni _loops ] n # [^-loops] \ ^ 

The first contribution is the product of the local bulk plaquette weights; the second contribution 
is the product of the local boundary plaquette weights; while the third contribution is the 
product of the non-local loop weights. Finally, the partition function is the sum over all possible 
configurations: 

Z = Y,P(G). (2) 

G 

Although we are working in a slightly different domain, the loop model with the boundary 
plaquettes just described was introduced by Dubail, Jacobsen and Saleur [26l|27]. We generalise 
this loop model in Section [3] and further in Appendix |A.l 
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Figure 3: A subdomain of the domain Q with the four types of loops shown: a bulk 
loop, of weight n, passing through bulk plaquettes only; a bulk loop, of weight n, 
passing through a /?2 boundary plaquette; a boundary loop, of weight 712, whose top 
endpoint lies in the bottom half of a /3s boundary plaquette; and a pair of boundary 
loops, of weight 774, whose top endpoints lie in top halves of /J3 boundary plaquettes. 
The acute angle of the rhombi is given by a. 

2.2 Parafermionic observable 

For the dilute 0(n) model with only bulk loops, Smirnov introduced a discrete, complex- valued 
observable defined as follows. 

One first introduces a defect consisting of an open, oriented loop segment starting at a base 
point, here set to be -00. A configuration 7 is then a collection of closed loops in addition to 
this defect. The observable is defined at each of the mid-edges of the lattice by 

F(z):= £ P(l)e- isW{l \ (3) 

where the sum is over configurations 7 for which the defect ends at z. P{^f) is the Boltzmann 
weight of the configuration 7, while W{^f) is the winding angle of the defect from the base point 
-00 to z. We define the initial winding angle to be along the positive real axis. The parameter 
s is referred to as the spin, giving rise to the term parafermionic observable. 

Smirnov introduced ^ with the aim of constructing a lattice observable that satisfies the 
discrete holomorphicity condition in the bulk: 

Y l F(z i )Az i = 0, (4) 

i 

where the sum runs over the edges of a plaquette or a cluster of plaquettes in a counterclock- 
wise manner. Ikhlef and Cardy [12] observed that Q precisely holds for the integrable 0(n) 
model, i.e. when the Boltzmann weights satisfy the Yang-Baxter equation. Here we recall their 
argument. 

The configurations appearing in the bulk dilute O(n) model arise from considering the four 
distinct ways a loop configuration external to a single plaquette can be fixed, see Fig. |U To 
illustrate, let us consider the case described in Fig.[5j and let P(G*) denote the total Boltzmann 
weight of the configuration external to the plaquette. The discrete contour integral arising from 
the local configurations in Fig. [5] is then given by 

£FCzj)Azj = P(G*)e- sin/2 [nu 2 + tf(nw 2 + wi) - - Cf4 (5) 

i 
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Figure 4: The bulk plaquettes of the dilute O(n) model: bulk connectivities. 









Figure 5: The possible internal connectivities for the external configuration indicated 
by the second diagram in Fig. |4j giving rise to the second equation in ([6]). 



From now on, we will discard overall irrelevant pre- factors such as P(G*). Thus, the following 
complex equations are obtained by requiring discrete holomorphicity Q, where the discrete 
contour integral follows from going around a single plaquette in a counterclockwise manner in 
each of the four cases of Fig. [IJ 

nv>i + C£ 2 ^ ~ i u 2 ~ C( w 2 + nwi) 
nu2 + (£,(nw2 +wi)- £ui - (^v 

2 

nv + C^ 2 ui - (C 2 ^i + £ w 2) ~ CC^2 

t + (£U 2 -V-CUl 

where 

£ = e is7r , C = e" ix , x = a(s-l). 
Assuming real Boltzmann weights, and parameterising the loop fugacity n by 



o, 
o, 
o, 
o, 



(6) 



(7) 



n = -2cos4A, 



(8) 



it follows from simple linear algebra that ^ only has a non-trivial solution if s = 3A/7T + 1. The 
Boltzmann weights that solve ^ precisely correspond to the well-known integrable weights 



t(x) = sinxsin(3A - x) + sin 2 A sin 3A, 

u\{x) = sin2Asin(3A - x), 

U2(x) = sin2Asinx, 

v(x) = sinxsin(3A - x), 

wi(x) = sin(2A - x) sin(3A - x), 

W2(x) = -sinxsin(A - x), 



(9) 



arising as the one-parameter family of solutions to the Yang-Baxter equation obtained in fT6^ H7j. 



3 Discrete holomorphicity with a boundary 



3.1 Bulk configurations 

The introduction of the (3$ plaquette results in a large number of additional bulk configurations 
as a loop segment may now pass through the boundary via these new boundary plaquettes. We 
now also have to modify the definition of the parafermionic observable to account for this and 
still be able to require that the observable satisfies the discrete holomorphicity condition in the 
bulk @. 






Figure 6: The bulk plaquettes of the dilute 0{n) model: boundary connectivities. 



Consider, for example, the case where the external loop configuration is attached to the 
boundary via a single /3s plaquette. In Fig. [6j we indicate this by placing a blob on the loop 
segment. In such cases, two things may happen. We may obtain a boundary loop weight 
n\ by closing the loop in the interior of the bulk plaquette. Alternatively, in case the defect 
becomes connected to the (3% plaquette, we say that it "passes through" the /3s plaquette before 
re-entering the bulk plaquette we are calculating the contour integral around. In such cases, we 
introduce an additional factor q\ {/q r ) to the parafermionic observable according to whether the 



defect has passed counterclockwise (clockwise) through the f3% plaquette, see (10) below. 






Figure 7: Defect passing through a /3s boundary plaquette before (and possibly again 
after) entering the bulk plaquette. 



A second set of additional configurations occurs if the defect segment has passed through a 
(3s plaquette before entering the bulk plaquette, see Fig. [7| To illustrate this, the configurations 
corresponding to the second connectivity in Fig. [7] are explicitly depicted in Fig. |HJ In this case, 
loops with weight ri2 may occur as well as a new type of boundary loop: a loop that passes from 
the top edge of a /3s plaquette to the top edge of another /3s plaquette, as shown on the left in 
Fig. |9j This type of loop is only possible because of the presence of the defect loop segment. 
We give such a loop the weight 713. Similarly, a boundary loop from the bottom edge of a (3s 
plaquette to the bottom edge of another /3s plaquette is also given the weight 77,3. 

For these reasons, we modify the definition of the observable to 

F(z):= £ Pbd y (7)e- W(7) (gi) e (g 1 f, e,e,£ + £€{0,l}, (10) 

7: — 00->Z 

where the Boltzmann weight Pbdy(7) is defined similarly to the expression for the weight P(j) 
of the configurations in ^ (see ([I])), but may also include contributions from the n3-type loops. 
The possible factor q\ or q 1 is associated to the defect if it passes through a /3s boundary 
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Figure 8: Detailed configurations corresponding to the second connectivity in Fig. ul 





Figure 9: The third type of boundary loop shown in red. It carries the weight n^. 



plaquette in a counterclockwise (e = 1, 6 = 0) or clockwise (e = 0, e = 1) manner, respectively. If 
the defect passes through more than one /3s plaquette, it is the orientation (counterclockwise or 
clockwise) of its last visit which determines the choice of factor q\ or q 1 . 
For the cases in Fig. [6j we thus obtain the following equations 

niui + q^^v - qiim - ((qiw 2 + niwi) = 0, 
n\u 2 + (£(niw 2 + - - ^lC?^ = 0, 

n\v + <ZiC£ V - {q x i 2 wi + q^ w 2 ) - qidm = 0. 

Writing 

q 1=p e 4iA+2iA i , ni = -2p cos 2 Ai , 

it is easy to verify that these are satisfied for arbitrary p with the integrable Boltzmann weights 
given in Q. The factor p is fixed below. 

For the cases in Fig. [7J we obtain the equations 

qiniui + C^q^sv - \q_\n 2 u 2 - (qi(n 2 w 2 + niwi) 
q\n\u 2 + C£ (91^1^2 + 9i^3^i ) " ^Zi^i - (Iqin 2 v 

2 

qiniv + C^q^ui - (^q^wi + £ gin 2 W2) - C,^qin 2 u 2 
Using the expressions Q for the Boltzmann weights, we then obtain a condition for n^: 

e 8Al g 1 (n 3 -gi) + gi(n 2 -gi) = 0. 



(11) 



(12) 



0, 
0, 
0. 



(13) 



Requiring that ns is real results in 



n 3 = -n 2 



sin 4 A 



sin(4A + 4Ai)' 



P = n 2 



sin2Ai 



sin(4A + 4Ai) 



(14) 



(15) 



or equivalently, 

n\ sin4Ai rt\ sin4Ai 

ri2 sin(4A + 4Ai)' n<$ sin4A ' 
from which it follows that the various loop weights are related by 

n 3 = n i + n 2 ~ nniri2. (17) 



We will comment on this relation in the discussion following (31). 



4 Discrete holomorphicity at the boundary 

Recall that, in the complex plane, we can write an integral of an analytic function F = F\ + 1F2 
over a contour C as the sum of two real- valued line integrals in M 2 , 

f c Fi,-f o F.di + if o F.dn, (18) 

where F - (Fi, -2*2 ), dl is the infinitesimal line element parallel to C, while dn is the infinitesimal 
outward normal to C. If C is a closed contour, holomorphicity requires both the real and 



imaginary parts in (18) to vanish. However, if we integrate along the boundary, C is not a 



closed contour. In this case, we require either the real part or the imaginary part of (18) to 
vanish. 

In the discrete setting, we therefore have two natural choices of boundary conditions for the 
parafermionic observable: 

Re(^ J P(^)A2 i )=0, (19) 

i 

if we require zero flux of F along the boundary, or 

Im(2F(^)A^) = 0, (20) 



if we require zero flux of F across the boundary. The sum in (19) and (20) is along the 
two edges of a single boundary plaquette, and counterclockwise orientation means along the 
edges in the downward direction. We will now show that requiring these boundary conditions 
separately (not simultaneously) results in integrable boundary weights, i.e. weights satisfying 
the (boundary Yang-Baxter) reflection equation [18] . Throughout, we shall assume that the 
boundary plaquette weights 02 and $3 are real. 



4.1 Boundary configurations 

The first set of boundary configurations to consider are the ones where the defect loop segment 
enters a boundary plaquette without having passed through a /3s plaquette prior to arriving at 
the boundary plaquette under consideration. When entering from below, as indicated to the 




Figure 10: First set of boundary configurations. 
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left in Fig. 10, the expression for the discrete contour integral is 

£ F( Zi )A Zi = r 1/2 (C/?i - - Cgi&) . 



(21) 



As we shall require that the observable satisfies the boundary condition (19) or (20), let us 
compute the real (R) and imaginary (/) parts of the above expression: 

Ri = sin(x - §A)/?i + sin(x + §A)/3 2 + psin(x - §A - 2Ai)/3 3 , (22) 
/i = -cos(x- §A)/3i +cos(x+ §A)/3 2 + pcos(x - §A-2Ai)/3 3 . (23) 

If the defect enters from above instead, we simply get minus the complex conjugate of (|2l|). 



Since our boundary condition, (19) or (20), amounts to setting either the real or the imaginary 



part to zero, the choice of entry edge is immaterial. This observation also applies to (24) and 



(|27j) below. 

As indicated in Fig. [TT] the second type of configurations corresponds to the defect entering 
the boundary plaquette after having already passed through a /3 3 boundary plaquette. For the 
case to the left in Fig. [TTJ we find 



i 

whose real and imaginary parts are given by 

R 2 = p[sin(x + |A + 2Ai)/3i +sin(x- |A-2Ai)(/3 2 + n 2 /3 3 )], 
h = P [ - cos(x + |A + 2Ai)/?i + cos(x - |A - 2Ai)(/3 2 + n 2 /3 3 )]. 



(24) 



(25) 
(26) 




Figure 11: Second set of boundary configurations. 

The final set of boundary configurations to consider involves those where the defect loop 
passes clockwise through a /3 3 boundary plaquette before entering the boundary plaquette under 



consideration in a counterclockwise manner (from below), or vice- versa. As follows from Fig. 12 



each of these cases includes the possibility of an n 3 -type boundary plaquette. Indeed, the discrete 
contour integral for the case to the left in Fig. [12] is given by 

(27) 



£ F(zi)A Zi = r 1/2 (C9ift " Cgi/?2 - Cqmzfo) , 



whose real and imaginary parts are given by 

R 3 = p[sm(x- ^A-2Ai)/3i +sin(x+ + 2Ai)/3 2 + n 3 sin(x - §A-2Ai)/3 3 ], (28) 
I 3 = /9[-cos(x-^A-2Ai)/3i+cos(x + ^A + 2Ai)/32 + n 3 cos(x-|A-2Ai)^3]. (29) 
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Setting Ri = R2 = R3 = yields three equations, one of which is linearly dependent on the 
other two, allowing us to solve for the ratios fii(x) / f3s(x) and ^(x)/ /3s(x). Likewise, we can 
solve for the same ratios in the case I\ = I2 = I3 = 0. This gives the solution sets 

Pi(x) = 713 [±cos(2x - A) - cos(4A + 4Ai)] , 

fo(x) = n 3 [±cosA-cos(2x-4A-4Ai)] , (30) 
(3s(x) = -2sin4A sin2x, 

where the + sign is for vanishing real parts, while the - sign is for vanishing imaginary parts. 
Written in terms of the boundary loop fugacities ni, 77,2 and 713 explicitly, the weights read 

Pi(x) = n\ + 712 cos4A ± ris cos(2x - A), 

fa (x) = n\ cos2x + ri2 cos(2x - 4A) ± 713 cos A, (31) 
/3s(x) = -2sin4A sin2x, 



where we recall the relation ttJ = n\ + n\ - nn\ri2 given in (17). 

We have verified explicitly that the weights ( |3T| ) precisely match the integrable boundary 
plaquette weights following from solving the (boundary Yang-Baxter) reflection equation for 



this model. From that perspective, 713 is merely a parameter defined by (17) since n3-type 
loops can not be formed by any loop segment but the defect. One can generalise the model, 
though, by introducing boundary plaquettes with a single loop segment, thereby allowing n<$- 
type loops to be formed. These new boundary plaquettes and the solution to the corresponding 
reflection equation are discussed in Appendix [ATj Whether a natural parafermionic observable, 
compatible with discrete holomorphicity, can be defined in this generalised 0(n) model remains 
to be seen. 



In the special case 



sin(4A + 4Ai) 
the boundary plaquette weights (|3T|) reduce to 



sm4Ai 

n i = — 7VT> ^2 = 1, (32) 



fii(x) = ^[cos(2x- A) =fcos(4A + 4Ai)], 

j3 2 ( x ) = ^[cosAtcos(2x-4A-4Ai)], (33) 
/3 3 (x) = ±sin(4A + 4Ai) sin2x, 



where we have introduced the sign- dependent rescaling 



ftM. ^^ aw. i- 1,2,3. (34) 

z sin 4A 
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In factorized form, these solutions are recognised as the ones appearing in the blobbed 0(n) 
model of Dubail, Jacobsen and Saleur [27J. Even though the relation n\ = n\ + n\ - nn\U2 is 



homogeneous in the boundary loop fugacities (ni, ri2 and 713), the specialisation (32) is not 
merely a rescaling (assuming n 2 ± 0) since /3s is not proportional to any of the loop fugacities. 
Our solutions are therefore more general than the ones in [27J. 

The diagonal or reflecting boundary conditions considered by Ikhlef [25] correspond to setting 



/3s = 0, in which case we should only consider the boundary configurations depicted in Fig. [10 
Imposing the boundary condition R\ - yields 

Pi(x) = sin(|A + x), (3 2 (x) = sin(|A - x), (35) 

while imposing the boundary condition I\ - yields 

0i(x) = cos(|A + x), (3 2 (x) = cos(|A - x). (36) 

As solutions to the reflection equation, these two solutions were originally obtained by Batchelor 
and Yung [23] . 

5 The loop model 

The loop model is a model of densely packed loops of two possible colours. In the bulk, 
a parafermionic observable has been defined by Ikhlef and Cardy [12]. After reviewing this 



definition in Section 5.1, we extend it to accommodate non-trivial boundary conditions. 

The set of bulk and boundary plaquettes are shown in Fig. [13} There are also three types 
of loops in the model. Loops that lie entirely within the bulk have fugacity n. As in the dilute 
O(n) model, loops that are attached to the boundary are of two types: those whose top end is 
attached to the top edge of a boundary plaquette have fugacity 774, while those whose top end 
is attached to the lower edge of a boundary plaquette have fugacity n 2 . The loop fugacities n 
and n\ are parameterised as 

n = -2cos4A, ri\ - -2pcos2Ai. (37) 
The parameterisation of 77,2, as well as how p is determined, is discussed below. 

5.1 Parafermionic observable 

We again introduce a defect in the model, this time consisting of a closed loop formed by gluing 



together two open loop segments of different colours, see Fig. 14 The value of the parafermionic 



observable F(z) is defined at the point z where the colour changes, 

F(z):= Y, P(~fX isW(l) . (38) 

The Boltzmann weight P(j) of the configuration 7 is defined as in the 0{n) model as the 
product of the plaquette and loop weights. The winding angle W(j) is the sum of the winding 
angles of the two loop segments of the defect, where the initial orientations of the segments are 
along the positive and negative imaginary axes, respectively. 

As before, the above definition has to be modified in the boundary model because of pos- 
sible boundary interactions of the defect loop. A configuration acquires a factor q\ (g x ) if the 
segment passes through the boundary in the counterclockwise (clockwise) sense. In certain con- 
figurations, the loop segment passes through the boundary more than once. Such configurations 
acquire a factor according to how the segment passes through the boundary the final time. 
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Figure 13: Bulk and boundary plaquettes of the loop model. The inclusion of 
the /3s and /?4 plaquettes gives rise to loop types not present in the bulk version of the 
model. 



Figure 14: The definition by Ikhlef and Cardy [12] of the parafermionic observable 
F(z) in the loop model. The open circle represents a point at -00, while the 
black disk is the point z at which the two colours meet. 



Owing to the presence of the defect loop, boundary loops can be formed between two upper 
edges or between two lower edges of a pair of boundary plaquettes. We assign such boundary 



loops the fugacity 713. As in the dilute 0(n) model (10), we thus define the parafermionic 
observable by 

F(z):= Y, ^bd y (7)e- W(7) (9ir(9i) e ", e,e,e + e e {0, 1}, (39) 

7: — 00— ^ 

where Pbdy(7) is l ne weight of the configuration 7 which may involve TT^-type loops. 

The parameters gi, 71 2 and 713 are fixed from requiring discrete holomorphicity in the bulk 



and that 713 be real. Here we summarise the results which are verified in Section [5T2j Thus, as 
in the dilute 0(n) model, gi, 712 and 713 are conveniently parameterised as 

^1 \(a\+9\i\ sin(4A + 4Ai) sin4A , . 

91 = " 9 \\ e l(4A+2Al) ? n 2 = -m \ U , n 3 = ni— p. 40 
2cos2Ai sm4Ai sm4Ai 

5.2 Bulk configurations 

The integrable bulk weights [28] are 
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Figure 15: The first set of bulk configurations in the C\ model. These configurations 
are the only ones appearing in the bulk model considered by Ikhlef and Cardy [12] . 




Figure 16: The first additional set of bulk configurations in the C\ model. In 
these cases, the defect passes through a fa plaquette, indicated with a blob, before 
re-entering the given bulk plaquette. Similar configurations with a fa plaquette follow 
by interchanging the two colours. 



u\{x) = sin 2 A sin(x - 6A) , 

U2(x) = - sin 2 A sin x, 

v(x) = -sinxsin(x - 6A), (41) 

wi(x) = sin(x - 2A) sin(x - 6A), 

W2(x) = -sinxsin(x - 4A). 



In the following, we use 



£ = e 2 ™, C = e" ix , z = a(2s-l), (42) 

3A _ 1 
2 2' 

First, we consider the configurations which also appear in the bulk model. These are indi- 



where s - 2 2 . 



cated in Fig. 15, and discrete holomorphicity gives the following equations 

nui + ££v - %U2 - n£wi - £w2 = 0, 

n%W2 + £wi + n(u2 - v - £(ui = 0, (43) 
nv + £C u i ~ £ w i ~ £ w 2 ~ C u 2 = 0. 
There are additional sets of bulk configurations, due to the possibility of the defect passing 



through the boundary via a fa or fa plaquette. In the first additional set, shown in Fig. [16 
the defect enters the given bulk plaquette before entering a fa or fa boundary plaquette. This 
gives the following three equations 

niui + q^Cv - q 1 lu 2 n 1 (w 1 - qi(w 2 = 0, 

ni^w 2 + qi£wi + niCu2 q x v - qi£(ui = 0, (44) 
^ + 9i£C^i "9i£wi "9^2 "9i(>2 = 0. 



In the second set of additional configurations, shown in Fig. [TTj, the defect passes through a fa 
or fa plaquette possibly twice: once before entering the given bulk plaquette and possibly again 
before re-entering. Requiring discrete holomorphicity imposes the conditions 

niqui + nsq^Cv - n 2 9i^2 rtiq^wi - n 2 qi(w2 = 0, 

nigi^ 2 + n 2 9i^i + m9iC^2-^39i^-^29i^C = 0, ( 45 ) 
niqiv + nsq^Cui - n^q^wi - n 2 9^2 ^29C^2 = 0. 
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Figure 17: The second additional set of bulk configurations in the model, here 
illustrated for a particular choice of colours. In these cases, the defect is blobbed before 
entering and possibly re-entering the bulk plaquette. 





4 




Figure 18: Another set of additional bulk configurations in the model, here 

illustrated for a particular choice of colours. In these cases, both coloured segments of 
the defect are blobbed; one of them before entering the bulk plaquette. 



It is straightforward to verify that both sets of equations are satisfied using the parameterisations 
of ni, ri2, ns and q\ in (37) and (40). This implicitly fixes the parameter p in (37). 

There are yet more sets of configurations to consider. One set arises from interchanging the 
colours of those already considered. However, due to the bulk weights being symmetric under 
interchanging colours, the equations from these configurations are readily satisfied. Other sets 
arise from configurations where a loop segment, that cannot close to form a loop, is blobbed, see 
Fig. [18) In all these cases, the ensuing equations are equivalent to the ones already considered. 



5.3 Boundary conditions and integrable weights 
5.3.1 Diagonal boundaries 

Diagonal boundary conditions correspond to /?3 = /?4 = 0, in which case we only need to consider 



the configurations in Fig. [T9J The discrete contour integral is given by 

£F(20A^ = C/?i-C&. (46) 

i 

The vanishing of the real or imaginary part gives, respectively, 

ft = ft, ft = -ft- (47) 




Figure 19: Boundary configurations in the case of diagonal boundary plaquettes. 
The black disk indicates the point where the two defect loop segments are joined. The 
defect loop segments originate from -oo, indicated here and in subsequent diagrams 
by an open circle. The blue (orange) loop segments are initially oriented along the 
positive (negative) imaginary axes. 
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Figure 20: The first set of configurations in the case of non-diagonal boundaries. 



Figure 21: The second and third set of configurations in the case of non-diagonal 
boundaries. 



5.3.2 Non-diagonal boundaries 

We now have four weights to solve for: /3i , /?2 ? /?3 and /J4, and five types of configurations to 
consider at the boundary. These are shown in Figs. 20, 21 and 22, In the first set (Fig. 20), 
the defect loop segments only pass through a single (3% or ^4 plaquette. The discrete contour 
integral is 

£ F(zi)Azi = -C/3i + C/?2 - qiCPs + 91C/84. (48) 

i 

The real and imaginary parts of the above expression are 

R\ = -cosx/3i + cosx/32 - pcos(x - 4A - 2Ai)/?3 + pcos(x - 4A - 2Ai)/?4, (49) 



Ii = sinx/3i + sinx/32 + psin(x - 4A - 2Ai)/?3 + psin(x - 4A - 2Ai)/?4, 



where 



P = " 



ni 



2cos2Ai 



(50) 



(51) 



In the second set, shown on the left in Fig. [21} the orange loop segment has already passed 
through a /3 s plaquette. The two loop segments form a closed loop with total winding angle 2tt. 
The discrete contour integral yields 



]T F(zi)Azi = -qiCPi + qiCfa ~ qiri2(/3 3 + q^iC/3^ 



(52) 



and its real and imaginary parts are 

i?2 = -pcos(x - 4A - 2Ai)/?i + pcos(x + 4A + 2Ai)/?2 - ri2pcos(x - 4A - 2Ai)/?3 + p 2 cosx/34,(53) 
7 2 = psin(x - 4A - 2Ai)/?i + psin(x + 4A + 2Ai)/3 2 + n 2 psin(x - 4A - 2Ai)/3 3 + p 2 sinx/3 4 . (54) 
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Figure 22: The fourth and fifth set of configurations in the case of non-diagonal 
boundaries. 



The third set, shown on the right in Fig. 21, is the same as the second, but the colours have 
been interchanged. This leads to the following discrete contour integral 

£ F(zi)Azi = giC/3i - qi(/3 2 + ?i9iC/% " ™2<?CAb (55) 

i 

whose real and imaginary parts are 

i?3 = pcos(x + 4A + 2Ai)/3i - pcos(x - 4A - 2Ai)/3 2 + p 2 cosx/3 3 - n 2 pcos(x - 4A - 2Ai)/?4 3 (56) 
7 3 = psin(x + 4A + 2Ai)/3i + psin(x - 4A - 2Ai)/3 2 + p 2 sinx/3 3 + n 2 psin(x - 4A - 2Ai)/?4. (57) 



In the final two sets of configurations, shown in Fig. [22J the defect has passed through the 
boundary before entering the boundary plaquette, but the two loop segments form a closed loop 
that crosses itself and has total winding angle 0. The resulting discrete contour integrals are 

£ F(zi)Azi = -q^fc + qiCfc - qlCPs + q^Cfa (58) 

i 

and 

£ F(zi) Azi = qiCp! - qi (p 2 + n&Cfo - ql(f3 4 . (59) 

1 

Taking real and imaginary components of each of the above gives 

i?4 = -pcos(x - 4A - 2Ai)/?i + pcos(x + 4A + 2Ai)/3 2 - p 2 cos(x - 4A - 4Ai)/3 3 

+ pn 3 cos(x - 4A - 2Ai)/?4, (60) 

J 4 = p sin(x - 4A - 2Ai + p sin(x + 4A + 2Ai )/3 2 + p 2 sin(x - 4A - 4Ai )/3 3 

+ pn 3 sin(x - 4A - 2Ai)/3 4 (61) 

and 

i? 5 = pcos(x + 4A + 2Ai)/3i - pcos(x - 4A - 2Ai)/3 2 + n 3 pcos(x - 4A - 2Ai)/3 3 

- p 2 cos(x - 8A - 4Ai)/3 4 , (62) 
Is = psin(x + 4A + 2Ai)/3i + psin(x - 4A - 2Ai)/3 2 + n 3 psin(x - 4A - 2Ai)/3 3 

+ p 2 sin(x - 8A - 4Ai)/3 4 . (63) 

Taking only the real expressions or only the imaginary expressions and setting each to zero 
(i?i = R2 = R3 = i?4 = i?5 = or I\ = I2 = / 3 = / 4 = I5 = 0), we obtain a linear system in 
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the unknown Boltzmann weights. One readily observes that only R2 and R3 are linearly 
independent and, similarly, I2 and ^3 are also linearly independent. These are, in fact, the 
only independent equations arising from the boundary. 

Now, considering the real equations and solving for the weights, we thus find 

fii{x) - ni cos(x - 4A + 4Ai), 

= fh(x), 

f3s(x) = 2sin4Aisinx, 

fa(x) = fo(x). 

Alternatively, setting the imaginary components to zero yields 

0i(x) = ni sm(x - 4A + 4Ai ) , 
fa{x) = -/3iO), 
(3s(x) = -2sin4Ai cosx, 
fa(x) = -fc(x). 

We have verified that both of these sets of solutions solve the corresponding reflection equation. 



6 Conclusion 

We have shown how to obtain integrable boundary weights in the dilute 0(n) and loop 
models by imposing simple boundary conditions on a discretely holomorphic parafermionic 
observable. Interestingly, in each model, two sets of solutions to the corresponding reflection 
equation are obtained, according to whether we require the real or imaginary part of the discrete 
contour integral of the observable along the boundary to vanish. 

In the case of non-diagonal boundaries, there are several types of loops, and we are forced 
to define new discretely holomorphic observables that incorporate additional parameters asso- 
ciated with boundary interactions. These parameters are fixed by requiring that (i) the discrete 
holomorphicity property in the bulk be preserved, and (ii) that a parameter which we interpret 
as the weight of a new type of loop be real. 

As is the case in the bulk, the discrete holomorphicity approach involves solving a simple 
linear system of equations and thus bypasses the complicated non-linear functional equations 
arising from the reflection equation. Once the solutions are obtained, it is straightforward to 
verify that they do indeed satisfy the reflection equation. This approach has allowed us to find, 
in a very direct way, new integrable boundary weights for both the dilute 0(n) and loop 
models. 

A natural extension of this work is to define discretely holomorphic observables in other 
geometries such as the infinite strip and the cylinder. In the case of the cylindrical geome- 
try, non-contractible loops appear, while in the infinite strip there are loops that touch both 
boundaries. 

We have only considered properties of discretely holomorphic observables, but it is also of 
interest whether expectation values of these observables can be evaluated in closed form. This 
was, for example, achieved in [29] for a particular observable in the Temperley-Lieb loop model 
on an inhomogeneous lattice and with open boundaries. 

Finally, a number of authors [HI [22j El] have found discretely holomorphic observables to 
be an important tool in proving the location of the critical points for self-avoiding walks. In 
the present work, focus has been on the integrable points of the models. However, since the 
integrable and critical points are expected to coincide in many cases, the observables defined in 
this paper could also be of use in rigorously proving the location of the models' critical points. 



17 



In [30], an off-critical extension was developed to explore what can be said rigorously about 
critical exponents. 
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A Reflection equations 



The Boltzmann weights of the bulk plaquettes are said to be integrable if they satisfy the Yang- 
Baxter equation. Similarly, the Boltzmann weights of the boundary plaquettes are integrable if 
they satisfy the reflection equation of the form 

Y, ^)W^(x + y)W;,(y)W^' (y-x) = £ W^ u (y-x)W;(y)W^,(x + y)w{ (x). (66) 
For our purposes, it is convenient to work with a diagrammatic representation of this equation, 



as shown in Fig. 23 The indices a, /3, 7 and 8 in (66) indicate the external connectivity, and 
we refer to them as the terminals. The other indices, which occur once on each side of the 
equation, are summed over. The weights are functions of a spectral parameter. 

To obtain an equation in the Boltzmann weights, one first specifies an external connectivity, 
and subsequently sums over all choices of Boltzmann weights such that the connectivity between 
adjacent plaquettes matches. Each term in the ensuing equation thus arises from a choice of 
bulk and boundary plaquettes that are consistent with the given external connectivity. 

From symmetry considerations, the only non-trivial equations arise from choices of external 
connectivity that are not invariant under interchanging the top and bottom terminals. For the 
O(n) loop model, there are thus five non-trivial functional equations arising from the reflection 



equation. The connectivities of the diagrams are shown in Fig. [24j As an example, Fig. [25 
shows the set of diagrams giving rise to 

03(y)v(x + y)fii{x)ui{x -y) + n 2 fi?>(y)ui(x + y)fo(x)v(x - y) 
+ h{y)ui{x + y)fi 3 (x)v(x -y) + fo(y)ui(x + y)/3 2 (x)v(x - y) 
= u 1 (x-y)f3 3 (x)v(x + y)f3 1 (y). (67) 

For the loop model, there are six non-trivial functional equations to consider, without 

assuming that the weights are symmetric under interchanging colours. The connectivities of 



the terminals are shown in Fig. 26 



A.l Generalised dilute 0(n) model 

One can extend the dilute 0(n) loop model described in Section [2] by including asymmetric 



boundary conditions, as indicated to the right in Fig. [27| In keeping with the notation for 
the other boundary plaquettes, we denote the weight of an asymmetric boundary plaquette 
by /?4. A crucial consequence of the introduction of these new plaquettes is the possibility of 
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Figure 23: A diagrammatic representation of the reflection equation. The spectral 
parameters associated with the plaquettes are shown. A functional equation in the 
Boltzmann weights is obtained for each choice of connectivity of the terminals, indi- 
cated here by open circles. 
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u 


o o 




o o 
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o o 


o o 



t 



Figure 24: The five connectivities leading to non-trivial functional equations for the 
dilute 0(n) model. The terminals are labelled by open circles. A black disk indicates 
that the associated terminal is connected to the boundary. 
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Figure 26: The six connectivities of the terminals leading to non-trivial functional 
equations from the reflection equation for the model. A blue (orange) disk indi- 
cates that the corresponding terminal is connected to a ^4 (/3s) boundary plaquette. 



forming boundary loops between two upper edges or between two lower edges of the boundary 
plaquettes. As we did for the parafermionic observable, we assign such loops the fugacity 713. It 
is stressed that a defect need not be present for such loops to appear in the generalised model 
discussed here. 

Assuming ^4 ± 0, we find that there is no solution to the corresponding reflection equation 
unless 

n x = n 2 = n 3 . 



In this case, 



/3i(x) - 2 cos Asin(|A + x) - k 2 n\ sin(^A + x) sin(^A - x) sin(|A - x), 

/3 2 (^) = sin(|A - x)[2cos A - k 2 n\ sin 2 (^A - x)], 

/3s(x) = -k 2 sin2Asin2xsin(^A - x), 

^4(x) = fcsin2Asin2x, 



(69) 



provides a one-parameter family of solutions, labelled by k. It follows, in particular, that there 
is no solution for which $4 but $3 - 0. 

We observe that, for k = 0, the solution (l69j) reduces to the solution (1351). We also observe 
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Figure 27: The four types of boundary plaquettes in the generalised dilute 0{n) 
loop model. The corresponding plaquette weights are denoted by /?2, Ps and /?4, 
respectively. 



that the solution approaches the solution ( |3T) ) (with ni = 712) in the limit 

lim - , p e {/3i,/3 2 ,/3 3 }; lim - = 0. (70) 

k^oo k 2 sin(2 A - x) k^oo k 2 sin(2 A - x) 
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